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EDGEWORTH EXPANSIONS FOR PROFILES OF LATTICE 
BRANCHING RANDOM WALKS 

RUDOLF GRUBEL AND ZAKHAR KABLUCHKO 


Abstract. Consider a branching random walk on Z in discrete time. Denote 
by Ln{k) the number of particles at site /c E Z at time n E Nq. By the profile 
of the branching random walk (at time n) we mean the function k i-A- Ln{k). 
We establish the following asymptotic expansion of Ln{k), as n ^ oo: 




e Fjjxnjk)) 

y/2Trip"{0)n ^ nJV2 



a.s., 


where r E No is arbitrary, (f{^) = log (^) cumulant gener¬ 

ating function of the intensity of the branching random walk and 


Xn{k) 


k — ip'{0)n 
V'¥’"(0)n 


The expansion is valid uniformly in A: E Z with probability 1 and the -Fj’s are 
polynomials whose random coefficients can be expressed through the deriva¬ 
tives of (p and the derivatives of the limit of the Biggins martingale at 0. 
Using exponential tilting, we also establish more general expansions covering 
the whole range of the branching random walk except its extreme values. As 
an application of this expansion for r = 0,1, 2 we recover in a unified way a 
number of known results and establish several new limit theorems. In particu¬ 
lar, we study the a.s. behavior of the individual occupation numbers Lnikn)-, 
where kn ^ depends on n in some regular way. We also prove a.s. limit 
theorems for the mode argmax^^^ Ln(k) and the height max/.^^ Ln{k) of the 
profile. The asymptotic behavior of these quantities depends on whether the 
drift parameter v?^(0) is integer, non-integer rational, or irrational. Applica¬ 
tions of our results to profiles of random trees including binary search trees 
and random recursive trees will be given in a separate paper. 


1. Introduction 

1.1. Statement of the problem. The branching random walk (BRW) is a model 
that combines random spatial motion of particles with branching; see m for a 
historical overview. In this paper, we restrict our attention to the BRW on the 
integer lattice Z. The model is defined as follows. At time 0 consider a single 
ancestor particle located at 0. At any time n S No every particle is replaced 
(independently of all other particles and of the past of the process) by a random 
finite cluster of descendant particles whose displacements w.r.t. the original particle 
are distributed according to some fixed point process C on Z. The number of 
particles in C may be random. The positions of the particles in ( need not be 


2010 Mathematics Subject Classification. Primary, 60G50; secondary, 60F05, 60J80, 60F10, 
60F15. 

Key words and phrases. Branching random walk, Edgeworth expansion, central limit theorem, 
profile. Biggins martingale, random analytic function, mod-</?-convergence, height, mode. 


1 








2 


RUDOLF GRUBEL AND ZAKHAR KABLUCHKO 


independent random variables (nor need they be independent of the number of 
particles). 

Denote by Nn the number of particles at time n G Nq. Then {Nni n G No} is 
a Galton-Watson branching process. We will always assume that this process is 
supercritical (meaning that m := EA^i >1) and with probability 1 never dies out 
(meaning that 7Vi 7 ^ 0 a.s.) The latter assumption could be removed, but then all 
results hold conditionally on non-extinction. Denote the positions of the particles 
in the BRW at time n by 

Our main object of interest is the occupation number Ln{k) defined as the number 
of particles located at site k G Z a.t time n € Nq: 

(1) L„{k) = #{1 < j < iV„ : Zj^n = k}. 

The random function L„ : Z —7 Nq will be referred to as the profile of the branching 
random walk at time n. The intensity of the branching random walk at time n is 
the measure on Z defined as the expectation of the profile: 

(2) ■dn{{k}) =V.Ln{k), kGZ. 

We will study the asymptotic shape of the profile L„ as n —>■ oo. More concretely, 
we will obtain an asymptotic expansion of L„ in powers of which is similar 

to the classical Chebyshev-Edgeworth-Cramer expansion for sums of independent 
identically distributed random variables. The latter will be recalled in Section [L^ 

It follows from the definition of the branching random walk that the intensity 
measure is the n-fold convolution of (so that, in particular, dniZ) = m"). 
Throughout the history of the BRW this fact has been used to relate the asymptotic 
properties of the random measures represented by L„, as n —)■ oo, to the analogous 
classical properties of convolutions of probability measures. In connection with the 
central limit theorem, this leads to the Harris conjecture P71 Chapter III, §16] which 
was proved (in various forms and for various models) in [371 UHl HI 1301 El [Ml [31 SO] ■ 
The work of Biggins maini refers to this analogy in connection with large deviation 
principles and local limit theorems. 

The present work adds one more example to this list by proving an asymptotic 
expansion of the profile. Since the asymptotic expansion contains many other clas¬ 
sical limit theorems as consequences, we will be able to recover many of the above 
mentioned results in a unified way (though under sub-optimal moment conditions). 
Our asymptotic expansion of the profile will be stated in Theorem |2.1| For 
comparison, an asymptotic expansion of the intensity measure dn will be given in 
Proposition [2.6| It turns out that the two expansions do not coincide: while the 
expansion of the intensity dn is deterministic, the expansion of the prohle con¬ 
tains random terms which can be expressed through the derivatives of a remarkable 
random analytic function given as the limit of the Biggins martingale; see below. 

Our motivation for investigating such asymptotic expansions was to develop a 
method to study binary search trees, random recursive trees, and some similar 
types of random trees that appear in the analysis of algorithms. An important 
characteristic of these trees is their profile, that is the random function k i—7 Ln{k) 
counting the number of nodes at a given level k, for a tree with n nodes. The 
probabilistic properties of these profiles have been much studied; see, e.g., [T8l 
EiiiniEij- Since random trees can be embedded into continuous-time branching 
random walks, see Chauvin et al. m and Biggins and Grey HU, it is possible to 
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translate our asymptotic expansion into the setting of random trees. In Sections |2.6| 
and |2.7[ we will answer the BRW analogues of several open questions on the profiles 
of random trees. For example, we will prove limit theorems on the height and the 
mode of the branching random walk. While we plan to give applications to random 
trees in a separate paper, we will point out connections to the existing literature 
on random trees here. 

The paper is organized as follows. In Section [TT2] we recall the classical asymptotic 
expansion for sums of i.i.d. random variables. In Section |1.3| we fix the notation 


and state our assumptions on the BRW. In Section 1.4 we briefly comment on 
our simulations. In Section]^ we state the asymptotic expansion and its numerous 
consequences. Proofs are given in Section 


1.2. Asymptotic expansion for sums of i.i.d. random variables. Let Zi, ^ 2 , • ■ • 

be independent identically distributed (i.i.d.) random variables with KZi = /r and 
Var Zi — > 0. Define the sequence of their partial sums: 

Sn = Zi -h . . . -h Z„. 


We assume that Z\ takes only integer values and that the lattice span of the distri¬ 
bution of Z\ is 1, that is there is no pair /i G {2, 3,...}, o S Z such that all possible 
values of Z\ are contained in the arithmetic progression a KL. 

Under these assumptions, the local limit theorem (see, e.g. Theorem 1 in [341 
Ch. VII, p. 187]) states that 

(3) lim i/nsup 

n-J-oo 

where 

(4) x^k) = kGZ. 

Under additional moment conditions, there is a complete asymptotic expan¬ 
sion of P[5'„ = fc] in powers of the classical Chebyshev-Edgeworth-Cramer 

expansion. Various versions of this expansion have been much studied; see the 
monographs by Petrov [33], Bhattacharya and Ranga Rao |3], Hall |32|. To state 
the expansion relevant to us, suppose additionally that E|Zi|’'+^ < oo for some 
r G Nq. The logarithm of the characteristic function of Z\ can then be written in 
the form 

logE[e*®^^] = ^ —I- o(|s|’'"''^), as s —>■ 0, 

i=i 

where the numbers Kj are the cumulants of Zx. Note that Ki = fj, and K 2 = <t^. The 
Chebyshev-Edgeworth-Cramer asymptotic expansion reads as follows, see Theorem 
13 in Petrov [331 Ch. VII, p. 205]: 


nSn = k]- 


\j2'Kn c 




= 0 , 


(5) 


r+1 

lim n 2 sup 

n-yoc 



kJ 


;(U 


Tin G 




i=o 


= 0 , 


where qj(x) is a degree 3j polynomial in a: = Xn{k) whose coefficients can be ex¬ 
pressed through the cumulants K 2 , ..., Kj+ 2 - The first three terms in the expansion 
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are given by 

2 

( 6 ) qo{x) = 1 , qiix) = ^Reaix), 92(2;) = ^^He4(x) + ;^He6(2:), 
where He„(a:) denotes the n-th “probabilist” Hermite polynomial: 

He„(x) = e5^^ 

The first few Hermite polynomials relevant to us are 

(7) Hei(x)=a;, He 2 (a;) = — 1, He 3 (a;) = a;^ — 3x, 

(8) He4(2:) =a:^-6a:2 + 3, He6(a;) = a;® - 15a:^ + 45x^ - 15. 

The aim of the present work is to obtain asymptotic expansions of this type for the 
profiles of branching random walks. 


1.3. Assumptions on the branching random walk. Consider a branching ran¬ 
dom walk on the integer lattice Z, as defined in Sect ion [lT| We will use the following 
standing assumptions. Recall that 411 denotes the intensity measure of the BRW 
at time n = 1. 


Assumption A: The branching random walk is non-degenerate, that is the set 
{fc S Z: 4?i({fc}) > 0} contains at least two elements. 


Recall that A„ denotes the number of particles in the BRW at time n. 
Assumption B: Ai ^ 0 a.s. and m ■= EAi e (l,oo). 

Assumption A simply excludes a trivial case. As we already mentioned in Sec¬ 
tion o the first part of Assumption B could be removed, but then all results hold 
conditionally on non-extinction. 

An important role will be played by the cumulant generating function p of the 
intensity i?i: 


(9) 


ip{l3) = logE 




e (-00, -1-00], 


Clearly, (/?(0) = logm > 0. 


/3 e K. 


Assumption C: The function (p is hnite in some open interval containing 0. 


Let Vip be the maximal open interval on which ip is finite. It follows that ip is 
strictly convex and infinitely differentiable on 'D^p. A crucial role in the study of 
the branching random walk is played by the Biggins martingale: 

(10) WnW) := /3 S V^. 

i—1 

By the martingale convergence theorem for non-negative martingales, we have for 
all /3 S T>p, 

(11) Wr^iP) ^ Woo(/3). 

n—¥(x> 

It is, however, possible that Woo(/3) = 0 a.s. The range of /3 where this does not 
happen was found by Biggins [5]. Denote by C the open interval on 
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which ip'{/3)P < 

(12) /3_ = inf{/3 e 

(13) /?+ = sup{/3 G : if’{13)13 < v?(/3)}- 


Clearly, 0 G (/3_,/3+), so that this interval is non-empty. The endpoints of the 
intervals and (/3_,/3+) are allowed to be infinite. We also need the following 
moment condition which supplements Assumption C. 

Assumption D; There is ap > 1 such that for every compact set K C (/?_,/?+), 

(14) sup E \Wf (/?)] < oo. 

IBeK 

Under the above assumptions, it follows from Biggins [3 Theorem A] that for all 
P G (/3_,/3+), 

Woo(/3) > 0 a.s., EWoo(/3) = l. 

It is a crucial observation due to Biggins and Uchiyama [31] that the 

martingale convergence 0 can be extended to a complex neighborhood of the 
interval (/3_, /3+). First of all, it is clear that the function (p is dehned as an analytic 
function of /3 in a sufhciently small open set U G C containing the interval V^p. It 
follows that for every n G N, the function Wn{P) is well-defined as an analytic 
function on U. It has been shown in [I] and |51| that under Assumption D there 
exist an open set V G C containing the interval (/3_,/3+) and a random analytic 
function Wao{P) on V such that 


(15) 


sup |IFoo(/3) - IU„(/?)| 
peK 


a.s. 



0 


for every compact set K G V. 

Assumption D is essential for our method of proof; see also Section 3.2 and the 
discussion in |8]. 

Finally, as in the classical Chebyshev-Edgeworth-Cramer expansion, we need to 
assume that the lattice width associated with the support of the intensity measure 
equals 1. 


Assumption E: There is no pair h G {2, 3,...}, a G Z such that i?! is concentrated 
on the arithmetic progression a + KL. 


Assumption E is not a restriction of generality because it can always be achieved 
by a suitable affine transformation; see Example |2.16 below. 


1.4. Simulations. At various places we illustrate our results by simulations. For 
Figuresand Figure]^ (top) these are based on the point process ^ given by 

P[C = M=P[C = <5-i + 5+i] = i, 

where 5z is the Dirac delta-measure at z G Z. In our simulations of the BRW 
we do not keep track of the individual locations as the number of particles grows 
exponentially. Instead we make use of the following rules to obtain L„: To(fc) = 
l{fc=o} and 


Ln+i{k) — {Ln{k) — Tn{k)) + Tn{k — I) + T„(fc -|- I), n G Nq, fc G Z, 
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where Tn(k) is the number of particles at time n and site k that have 2 descendants 
at time n + 1 (located at sites fc + 1 and k — 1). Given Ln{k), fc G Z, the random 
variables Tn{k), fc G Z, are conditionally independent and T„(fc) ^ Bin(L„(fc), 1/2). 

For Figurej^ (middle, bottom) and Figuresthe simulations are based on the 
point process C with 


P[C = <5o]=P, F[C = 2S+1 + S-i] = 1 - p 
with p = 1/2 (rational case) and p = 3(7r — 3) (irrational case). 


2. Results 

2.1. Asymptotic expansion of the profile. Consider a branching random walk 
on Z which satisfies Assumptions A-E. Recall that A„(fc) denotes the number of 
particles of the branching random walk that are located at site fc G Z at time 
n G Nq. Take some /3 G (/3_, /?+). Our limit theorems will be stated in terms of the 
“tilted” prohle 

fcHAe'^'=-‘^(^)"L„(fc), fcGZ. 

As known from large deviations theory, the tilting operation allows to better ac¬ 
cess the properties of the profile for fc « ip'{j3)n. Define the corresponding tilted 
cumulant 

(16) KjiP) = 

as the j-th derivative of ip at /3, j G N. In particular, we need the notation 

(17) p(/3) := (p'(/3) = ki(/ 3), cr^(^) := p"(/3) = N2(^) > 0. 

Introduce the “standardized coordinate” 

(18) X = Xnik) := ^ ’ fcGZ. 

a{/3)^/n 

Note that for ease of reading we omit the dependence on /3 in Xn{k). 

Our first theorem gives an asymptotic expansion of the tilted occupation number 
e^^“‘^(^)”'L„(fc) in powers of 

Theorem 2.1. Consider a branching random walk satisfying Assumptions A-E. 
Fix r G No and a compact interval K C (/3_,/3+). Then we have 


(19) 


u + 1 

n 2 sup sup 
peK fcez 


e'^'=-^(^)”L„(fc) 


e-k-liC ” 
a{l5)V2Trn ^ 



{xn{k);l3) 



n—¥oo 


where Fj{x;P) is a degree 3j polynomial in x. The coefficients of Fj(x; (3) are 
random and can be expressed through 

k2(/3), ..., Kj+2il3) and W^{l3),W^{/3),W^\l3). 

Stating the complete formula for Fj{x;P) requires introducing complicated no¬ 
tation and is therefore postponed to the proof in Section [373| Here, we provide only 
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the first three terms of the expansion: 

Fo(a:;/ 3 ) = t^oo(/ 3 ), 


( 20 ) 

( 21 ) 

( 22 ) 




F2{x-,j3) = WM 


6(t3(/3) 

^4(/3) 
24ct4(/3) 


He4(a;) 


4m 


He6(x) 




w'M 


720-6 (;3) 
1 


2a2(/3) 


He2(a:). 


Remark 2.2. We will see in the proof of Theorem 2.1 Equation 


that 


Fj{x-,P) = 


3 

E- 

771=0 




Qm,j [x'l 4) 1 


where Qo,j(x; /3),..., Qjj{x] /3) are deterministic polynomials in x with coefficients 
depending on /3. We will also see that Qoj{x;P) (the coefficient of Woo(,d)) is the 
same as qj{x) in ([^ but with Kj replaced by Kj{/3). 


Remark 2.3. Theorem 2.1 remains valid if in the formula for Fj we replace the 
derivatives of Woo by the corresponding derivatives of Wn- The reason is that 
w.p. 1, Wn converges to Woo exponentially fast together with all its derivatives 


(see Lemma 3.3), while the error term in Theorem 2.1 is of polynomial order only. 
This fact is used in our simulations where we replace Woo by Wn- It was shown 
in [35], see also [53], that in a suitable range of (3 the asymptotic distribution of 
the appropriately normalized difference Wn{P) — Woo{P) is a mixture of centered 
normals. A functional limit theorem for this difference was obtained in 1251. 


Remark 2.4. For the branching Brownian motion, an expansion similar to (19) 
(with /3 = 0) was obtained by Revesz et al. |3S]. In a general BRW, the transition 
mechanism is not Gaussian, so that methods specific to the Gaussian setting cannot 
be used. 

In the rest of Section we state numerous consequences of Theorem |2.1| It 
allows us to recover many known results in a unified way and to answer a number 
of open questions on the individual occupation numbers L„(fc), as well as the height 
and the mode of the BRW profile. 

2.2. Local and global central limit theorems. Taking only the first term in 
the expansion given in Theorem |2.1| (meaning that r = 0) we obtain 


(23) 


/n sup 

k^Z 




\/27r(p"(/?)n 


exp 


{k - 4mn) 

2(/?"(/3)n 


0 . 


The most interesting case is /3 = 0 when there is no tilting and we obtain the 
following local limit theorem for the BRW: 


(24) 


/nsup 

kGZ 


Lnik) Woo(O) 


427r<f"{0)n 


exp < - 


(fc — (p'{0)ny 
2ip"{0)n 


0 . 


Roughly speaking, (24) says that the profile k i—)■ Ln{k) has an approximately 


Gaussian shape with mean (^'(0)n, standard deviation 4 :/j"(0)n, and the total 
mass of the profile is iV„ ~ Woo(0)m"’; see the left part of Figure In Figures 
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and 1^ we use vertical bars for the discrete profiles and continuous lines for the 
approximating functions. 




Figure 1 . Left: A realization of the profile k i-A Ln(k) of a branch¬ 
ing random walk (blue) at time n = 50 and the first, Gaussian term 
(see ( 24 l) in the asymptotic expansion (red). Right: The same profile 
(blue) and the Gaussian term shifted by Wn{0)/Wn(0) (red). 


Let us now look at the next term in the expansion. Theorem |2 .1 1 with r = 1 and 
13 = 0 yields that 


(25) 


L„(fc) 

LFoo(0)m" tT(0)\/27rn 
cr(0)V27rn y/n V6(t3(0) 


He3(a:„(fc)) 


W^jO) Xnjk) 

fFoo(O) a(0) 


a.s., 


where the o-term is uniform in A: € Z. There are two correction terms on the 


right-hand side of (251. The term involving Hea is the “shape correction” to the 


Gaussian profile. The same term appears in the expansion of the expected profile 


EL„(A;); see Section 2.3 The term involving W!^{0)/Waoi0) can be thought of as 
a random “location correction”. Indeed, this term says that in order to obtain 
a better approximation to the BRW profile we have to take a Gaussian profile 
centered at '^'{0)n -\- W^{0)/Woo{0) rather than at <p'(0)n; see the left and right 
parts of Figure 

Note that the random variable W'^{0) appearing above is the a.s. limit of the 
martingale 

1 


= ZTf - V3'(0)n). 


Since the total number of particles at time n is Nn = IF„(0)77i"', we can view 
W'j^{0)/Wn{0) as an estimate for the “shift” of the profile w.r.t. its “expectation” 
ip’{0)n. This explains the appearance of W’^[0)/Wao{0) as a “location correction” 
in a quite natural way. Similarly, the variable (0) which will appear frequently 
below is the limit of the martingale 


N, 
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There is also a global central limit theorem for the BRW, originally known as 
the Harris conjecture HZl Chapter III, §16]. It states that for all a: G K, 


(26) — 
mr 


Zt,n - (p'iO)n 


< X 


^ IHoo(O) 

n—>-oo 


a/Stt J-c 


dy. 


Various forms of (261 and (241 have been obtained in [371 HSl dl 1301 El 1301 IB SB 123] ■ 
In fact, we can obtain a full asymptotic expansion in (26). To this end, one takes 


sums in Theorem |2.1| and uses the Euler-MacLaurin formula to approximate sums 
by integrals. We will record here only the first non-trivial term of this expansion. 

Proposition 2.5. Consider a branching random walk satisfying Assumptions A-E. 
Then, 


(27) 


1 


fCoo(0)r 




Ak) 


a/StT J-c 


dz = 


h— — oo 

/I 


_ «3(0) 
a{0)y/2Trn 6(T^(0) 

where the o-term is uniform over k € Z. 


bEoo(O) 


A non-lattice version of Proposition 2.5 has been obtained recently in 


In 


fact, similar results hold for many models other than the BRW; see [29] . 

2.3. Comparing the profile and the expected profile. It is interesting to com¬ 
pare the expansion of T„(fc) stated in Theorem 2.1 with the expansion of EL„(fc). 
Since the intensity measure of the branching random walk at time n is just the 
n-fold convolution of i?i, we can apply the classical expansion <§ to the expected 
profile EL„(fc) = d„({fc}). The proof of the following proposition is standard and 
will be given in Section |3.2[ Formally, it can be obtained from Theorem |2.1| by 
taking the expectation and noting that EIToo(/3) = 1 for all f3 G (/3_,/3+) which 
implies that all higher derivatives of Woo{J) have zero expectation. 

Proposition 2.6. Consider a branching random walk satisfying Assumptions A, 
C, E. Fix r G No and a compact interval K C Pip. Then we have 

— -x^ {k) ^ -| 

(28) + , 

where qj{x; /3) is the same as qj{x) in ([^ but with Kj replaced by Kj(/3). The o-term 
is uniform in /3 € K and fc G Z. 


We can considerably simplify the expansion given in Theorem |2.1| if instead of 
Ln{k) we consider Ln{k) — Woo(/3)EL„(A:). This normalization removes all terms 
involving lToo(/3) and we obtain 

Corollary 2.7. Consider a branching random walk satisfying Assumptions A-E. 
Fix r G No and a compact interval K C {l3-,/3+). Then, 

(29) e^'=-‘"(«"(L„(fc) - Woo(/3) EL„(fc)) 
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where F°{x]^) is the same as Fj{x;^) but without the term involving Woo{^), that 
is 


(30) 


FVx;l3) := Fj{x;(3) - Woo{l3)q]{.x\ (3). 


The o-term in (29) is uniform in (3 G K and k £ 7 a. 


Note that the summation in (29) starts with j = 1. We will need the first two 
terms of the expansion, c.f. (21) and (22): 


(31) Ff{x;/3) = ^^WUI3), 

(32) F°{x-J) = (/3)He2(a:). 


2.4. Uniform expansions. For every point k G 7 Theorem |2.1| yields a family of 
asymptotic expansions of Ln{k) parametrized by /3. It is natural to choose /? = /3(^) 
as the solution to ^ because then Xn{k) = 0. The information function 

I is defined by 

With the above choice of f3, Theorem |2.1| and Proposition |2.6| yield the following 
result. 


Corollary 2.8. Consider a branching random walk satisfying Assumptions A-E. 
Fix r G Nq. Then we have 


(33) e^F^)L„{k) = 


cr(/3(|))V^ 


1=0 


a.s., 


where the o-term is uniform over all k G 7 for which /3{^) exists and stays in a 
fixed compact set K C (/3_,/3+). Similarly, 


(34) 


^^(S)EL„(fc) = 






where the o-term is uniform over all k G 7 such that /?(-) exists and stays in a 
fixed compact set K C V^p. 


The expansion of EL„(fc) is well known because it is a classical formula for the 
precise large deviations of sums of i.i.d. random variables; see [T2l[2l|33]. Note that 


the terms with odd j in (33) and (34) vanish; see Remark 3.6 below. 


Remark 2.9. By taking a linear combination of (33) and (34) we obtain an expansion 
of 




L, 


.(fc)-lUoo(/3(^))EL„(fc) 


which looks exactly as (33) but with Fj replaced by F°. 


Taking r = 0 in the above theorem yields the following known result pm ; see 
Figurej^ In the context of random trees, results similar to Corollary |2.10| have been 
obtained by Chauvin et al. m. Drmota et al. [ID], Chauvin et al. [H Theorem 3.1] 
and Sulzbach [55] . 



























EDGEWORTH EXPANSIONS FOR PROFILES OF LATTICE BRW’S 


11 



Figure 2. A realization of the log-profile A: i—^ log Ln ( k ) of a branch¬ 
ing random walk (blue), the information fnnction k i—>■ (black), 

and the approximation obtained by taking the logarithm in (351 (red) 
at time n = 50. The realization is the same as in Fignrej^ 


Corollary 2.10. Let K be a compact subset of the interval (/3_,/3+) (for (HI;, 
resp. T>^p (for (36);. Then, 

(^(f))s"'<i>L„(l)-H„(;3(f)) ^^0. 


(35) sup 

feGZ: P{^)GK 


(36) 


sup 


(/3(0)e"^(*)EL„(A)-l 


0 . 


If fc = A;„ is an integer sequence such that lini„_>.oo kn/n = ^'{(3) with some 
/3 g (/3_,/3+) (for ( [^ ) ov (3 € V^p (for ([38|), then —>■ (3 and we obtain that 


(37) 

(38) 


n—^oo 

V2^(^"(/3)ne’^^(‘^)EL„(fc„) 1. 


Equation 


of Corollary 2.8 holds uniformly in fc G Z as long as - stays in 


the interval (i^'(/3_), i^'(/3+)) and remains bounded away from its ends. Here, both 
derivatives are understood as the corresponding one-sided limits and are allowed 
to be infinite. One may ask whether the whole range of the BRW is covered or 
just some part of it. A basic result due to Biggins [3], see also nni, describes the 
asymptotic range of the branching random walk: 


(39) 
where 

(40) 


1 a.s „ 

- max Zin — l-L, 
n i=l,...,Af„ n—>oo 




1 

— mm 


p 

^2,n ' ^ — 5 

n—¥oo 


T- = sup < (/5'(0). 
/3<o P 


If we additionally assume that ip is finite everywhere on M (rather than on some 
interval) or that the support of the intensity measure is finite, then it is easy 


to prove (see, e.g., ^ 
tion (p3| of Corollary 


) that p'{f3-) = r_ and p'{l3-^-) = r+. In this case, Equa- 


2.8 covers the whole range of the “central order statistics” 


of the BRW. The “extremal” and “intermediate” order statistics at r±n=Fo(’i-) 
are not covered. However, in some (rather exotic) examples it is possible that 
(p'{f3-),pfl3+)) is a strict subinterval of (r_,r+). 
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Example 2.11 (see [33]). Consider a BRW for which -did/c}) = ce~^/k^ for fc G N 
and i?i({fc}) = 0 otherwise. Here, c > 0 is a parameter. Then (p(/3) is finite for 
/3 < 1 and moreover, the left derivative Lp'{!—) := lim^^i is also finite, whereas 
for /3 > 1 one has <p(/3) = +oo. 



If c is sufficiently large, then the BRW is supercritical and we have 

This means that /3+ = 1, whereas r_|_ = (/^(l). Hence r+ is strictly larger than 

2.5. Continuous-time branching random walks. All results of this paper ap¬ 
ply to continuous-time branching random walks on Z which are defined as follows. 
At time 0 one particle appears at position 0. After an exponential time with param¬ 
eter A > 0, the particle disappears and at the same moment of time it is replaced 
by a random cluster of particles whose displacements w.r.t. the original particle 
are distributed according to some fixed point process ^ on Z. The new-born parti¬ 
cles behave in the same way as the original particle. All the random mechanisms 
involved are assumed to be independent. Denote the number of particles at time 
t > 0 hy Nt and note that {Nt : t > 0} is a branching (Markov) process in con¬ 
tinuous time; see m Chapter V]. Note that the law of the continuous-time BRW 
is uniquely determined by the following two parameters: the intensity A and the 
law of the point process C. The occupation number Lt{k) is defined as the num¬ 
ber of particles located at site /c G Z at time t > 0. If we restrict the time t to 
integer values only, we obtain a discrete-time BRW called the “discrete skeleton” 


of the original continuous-time BRW. If Assumptions A-E of Section 1.3 hold for 
the discrete skeleton, then all the results of the present paper can be translated 
in an evident way to the continuous-time setting by replacing n G Ng by t > 0. 
Note, however, that one has to be careful whenever the arithmetic properties of 
ip'{0) are involved; see Sections 2.6 2.7 The proofs require only straightforward 
modifications. 


2.6. Strong limit theorems for the occupation numbers L„(A:„). Recall that 
L„(fc) denotes the number of particles located at time n at fc G Z. Let us take an 
integer sequence fc = which behaves in some regular way. We ask whether 
the random variables L„(fc„) have a non-degenerate a.s. limit, after an appropri¬ 
ate affine normalization. The next proposition is known and follows immediately 
from (j^l and ([3^. 


Proposition 2.12. Consider a branching random walk satisfying Assumptions A- 
E. Let kn be an integer sequence such that kn ~ (p'{/3)n for some /3 G (/3_,/3+). 
Then we have 


An(^n) 

EL„(fc„) 


n—>-oo 




( 41 ) 
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Next we ask whether we can obtain more refined limit theorems for the “cen¬ 
tered” variables 


(42) (fc„) := L„(A:„) - W^{p) EL„(A:„). 

This question is especially natural if /3 = 0 and any particle in the BRW generates 
the same number m £ {2,3,...} of descendants. Then Woo(O) = 1 and hence 
the limit random variable provided by Proposition |2.12'| is a.s. constant. The same 
phenomenon occurs in the setting of random trees, where the natural analogue of 
bkoo(O) is equal to 1. 

We consider an integer sequence which, for some f3 £ (,0_,/?+), is represented 
in the form = ip'{j3)n + Cn- The result will depend on the asymptotic behavior 
of Cn- Recall that = p''{l3)- 


Theorem 2.13. Consider a branching random walk satisfying Assumptions A-E. 
Let kn be an integer sequence such that 

kn = p'iP) n + Cn 

for some /3 £ (/3_,/3+) and some Cn = 0{y/n). 

(a) If Cn = aa(fj)y/n -L o{y/n) for some a £ K, then 

(43) ne^'="-‘^(«"L“(fc„) ^ ^ 

n-s-oo y/2Tra^{fi) 

(b) If lim„_>.oo Cn = -too, but Cn = o{\/n), then 

(44) n3/2c-ie^'=-^(«"L°(fc„) ^ W'M- 

1^^°° V 2710-3 (/3) 

(c) If Cn is bounded, then 

(45) (fc„) - R(c„) ^ 0 

n—>-oo 

where 

(46) R{c)-.= -^ - ( (c+ - x^"(^)l . 

^ ^ v^a3(/3) V V 2a2(/3)y 2 

In the setting of random trees, Fuchs et al. [52] showed that the analogue of the 
occupation number L„(fc„), centered by its expectation and normalized by standard 
deviation, does not have a limit distribution fc„ = p'{fi)n + 0(1). Part (c) of the 
above theorem clarifies the structure of the set of a.s. subsequential limits in the 
BRW setting. 


Remark 2.14. 
quence 


In the situation of Theorem 2.13(c), limsup and liminf of the se- 


(47) 


^3/2g/3fe„-(^(/3)n 


L°n{kn) 


are a.s. finite (but not necessarily equal to each other). Whether or not the se¬ 
quence (47) has an a.s. limit depends on certain arithmetic issues; see Figurej^ We 
consider the following cases. 


Case 1. If the “drift” p'{l3) is integer, it is natural to take kn = p'{fi)n + a for 
some constant a £ Z. In this case, c„ = a which implies that R{cn) converges a.s. 
to R{a) and so does the sequence (47). This case is shown on top of Figure 


Cases 2 and 3. If p'{l3) is non-integer we cannot choose as in Case 1. Instead, it is 
natural to take = \}p'{j5)n\ +a, where a G Z, which means that c„ = a—{p'{P)n}. 
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Figure 3. A realization of the sequence L'^{kn) as a function 

of n = 1 ,..., 200 for = Li^'(0)nj + a with a = —1, 0,1 (red diamonds, 
green circles, blue triangles, respectively). The four black horizontal 
lines show the values R{—2),R{—1),R{0),R{1). Top: ¥’^(0) is integer, 
middle: V5^(0) = | is rational, bottom: ^5^(0) is irrational. 


Here, [-J denotes the floor function and {•} denotes the fractional part. The limit 
of c„ (and hence the limit of i?(c„)) does not exist. Instead, we can describe the 
set of all a.s. subsequential limits of (47). 


Case 2. For non-integer rational ‘f'iP) = | with coprime p and q the set of a.s. 


subsequential limits of (47) is finite and given by 


(48) 


: j = 0, ...,g- l| . 
In fact, over the subsequence u; = j + Iq, I € No, we have 


n3/2e^'="-‘^(«”L°(fc„) ^ r( a- 

n—>-OQ \ 
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This case (with g = 4) is shown in the middle of Figure]^ 

Case 3. For irrational (/?'(/?) the set of a.s. subsequential limits of (47) can be 
continuously parametrized by the interval [0,1]: 

(49) {i7(a-z) : ze [0,1]}- 

This case is shown on the bottom of Figure 

Remark 2.15. The above Cases 1 and 2 apply to BRW in discrete time only. In the 


case of continuous time t > 0, see Section 2.5 it is natural to take kt = [ip'{l3)t\ +a 
with a G Z. We have Ct = a — {gj'(/3)t} and there are two cases. 


Case 1. If = 0, then Ct = a and the sequence (47) has an a.s. limit R{a). 


Case 2. If (p'{/3) ^ 0, then the set of limit points of Cj is the interval [a — 1, a] and 
the set of a.s. subsequential limits of (47) is given by (49). 


Example 2.16. Consider a simple symmetric branching random walk on Z. That 
is, if at time n some particle is located at fc G Z, then at time n + 1 it moves to 
one of the sites fc — 1 or fc + 1 with probability 1/2 and generates there a random 
number of offspring. The distribution of the number of offspring is assumed to have 
finite mean m > 1 and finite p-th moment, for some p > 1. Note that <p^(0) = 0 by 
symmetry. Denote by zi ^, ■ ■ ■, Zn,Nru positions of particles at time n. Note that 
although Assumption E fails to hold for the simple symmetric BRW because the 
intensity measure is concentrated on {—1,1} C 2Z + 1, the transforme d BR W 
+n), 1 < i < Nn, satisfies Assumptions A-E. Applying Theorem 2.17(c) 


to the new BRW and switching back to the simple symmetric BRW, we recover a 
result obtained by Chen m Theorem 2.1] under a second moment condition on 
A^i. Parts (a) and (b) of Theorem 2.17 describe further possible asymptotic regimes 


that are not covered by m- Our Proposition |2.5| generalizes Theorem 2.2 of [T3] 
to more general BRW’s. Very recently, Gao and Liu elaborated on the result 
of Chen |15j by relaxing his second moment condition on Ni, allowing for general 
displacements distributions and by considering branching random walks in random 
environment. Note that Gao and Liu [23j impose a strong non-lattice assumption 
on their BRW’s and so their results do not cover the lattice case considered here. 
In particular, Gao and Liu [53] obtained non-lattice versions of Theorem 2.17 (c) 


and Proposition 2.5 


2.7. Height and mode of the branching random walk. Define the height M„ 
and the mode Un of the branching random walk profile at time n by 


(50) 


= maxL„(fc), = argmaxL„(fc). 


In the context of profiles of random trees, similar quantities (called width and mode) 
were studied by Chauvin et al. [T3|, Katona |3T], Drmota and Hwang [T9| and De- 
vroye and Hwang m- Here, we will consider the more general setting of branching 
random walks; applications to random trees are postponed to a separate paper. It 
is clear from the approximate Gaussianity of the profile, c.f. (24), that the mode 


should be near (p'{0)n, while the height should be approximately to"/ {y/2Trn cr(0)). 
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Theorem 2.17. Consider a branching random walk satisfying Assumptions A-E. 

(a) There is an a.s. finite random variable N such that for n > N, the mode Un 

is equal to one of the numbers [m*J or where 

^^(0) ^3(0) 

ITooW 2^2 (0) 

and [-J, [•] denote the floor and the ceiling functions, respectively. 

(b) The height Mn satisfies 


(51) 


< = T'i0)n- 


r(0) M„ 


= 1 - 


1 


(5^ ■ + (logM^oo)"(0)) + o 


Woo(0)to" ^ 2cr2(0)n \^6cr4(0) 4cr2(0) 

a.s., where 0n := min^gz |m* — k\ is the distance between u* and the closest integer. 
The second part of the above theorem obviously implies that 


(52) 


lim cr(0) Mn 


m" 


W'oo(O). 


In the context of random trees, similar results were obtained in [13 HU [HUT]. It 
has been asked by Drmota and Hwang [191 Section 5] whether M„, after appropriate 
centering and scaling, converges in distribution, and whether the limit distribution 
is (in our notation) the distribution of the logarithmic derivative Wf,{Q)/Waoifii)- 
Note that the analogue of Woo (0) equals 1 in the setting of random trees, so that ( |^ 
gives an a.s. constant, degenerate limit. Theorem 2.17 (b) answers these questions 
in the context of the BRW. In fact, we even have a limit theorem on the a.s. behavior 
of Mn. In our result, the second (rather than the first) logarithmic derivative of 
IToo at 0 appears. 


Remark 2.18. The first part of Theorem [TTT] does not say which of the numbers, 
[Un\ or [w*], is the mode. From the proof (which will be given in Section 3.6) 
it will be clear that for every given £ > 0 we can choose the random variable N 
such that for every n > N satisfying min^gz — k — \ \ > e the mode is equal 
to Un rounded to the nearest integer. In the case when is sufhciently close to a 
half-integer fc -I- k G Z, further terms of the asymptotic expansion are needed to 
tell whether u„ is [u* J or [m*] . It is natural to conjecture that the random variable 
^oo(0)/W^oo(0) is absolutely continuous. We would then have the following cases. 
Case 1. If v?^(0) is integer, then either Un = for all n > IV or Un = for 
all n > N. Indeed, the fractional part of m* is some constant which is not equal to 
1/2 w.p. 1. 

Case 2. If Vj'(O) is non-integer rational, then for n > N, Un chooses its value among 
[Un\ and I’m*] periodically; see the top of Figure]^ which shows the sequences [m*J, 
[m*] (both black circles) and m„ (red dots). 

Case 3. If <p'{0) is irrational, then there is no periodicity and both values [m* J and 
[m*] are attained with asymptotic density see the bottom of Figure]^ 

Remark 2.19. Since |0„| < 1/2, limsup and liminf of the sequence 

\/2'!Tn (t(0) Mn\ 


(53) 


Mn := 2a^{0)n 1 - 


IFoo(0)m" 
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Figure 4. A realization of the mode Un (red dots) as a function of 
time n = 1,... ,50. Black circles show and \Un~\- Top: The drift 
1 / 5 '(0) = I is rational, periodic behavior. Bottom: The drift is 

irrational, no periodicity. 


are a.s. finite. Whether or not the sequence (53) has an a.s. limit, depends on the 
arithmetic properties of Let 


(54) 


Q{0) = 


4(0) 

6cr4(0) 


K4(0) 

4cr2(0) 


+ 02 + (logWoo)"(O). 


Case 1. If :/9'(0) is integer, then On — 9 \s constant and the a.s. limit of (53) exists 
and equals Q{9). 

Case 2. If :/3^(0) is non-integer but rational, we have finitely many a.s. subsequential 
limits of On and of (53). Hence the normalized height (53) has finitely many a.s. 


subsequential limits inside the interval [Q(0),Q(5)]; see the top of Figure]^ Note 
that these limits are not equally spaced because of the quadratic term 0"^ in ( [54| . 
Case 3. For irrational the set of subsequential limits of On is the whole 

interval [0,^]. The set a.s. subsequential limits of (53) is therefore the interval 
{Q{9 ): 0 S [0, i]}; see the bottom of Figure]^ It is an interval of length 1/4. On 
Figurej^one sees that there are much more points close to the lower limit <3(0) than 
to the upper limit Q{\)- To explain this, note that the sequence is uniformly 
distributed on [0, |] by Weyl’s equidistribution theorem (stating that the sequence 
of fractional parts {(p'(O)n}, n S N, is uniformly distributed on [0,1]), hence the 
asymptotic density of 9“^ equals da;/\/^ on [0, i]. 


Remark 2.20. For branching random walks in continuous time t > 0, see Section [2.5[ 
we have the following two cases. 

Case 1. If ^^'(0) = 0, then 


0 := On = min 
feGZ 


wLio) 

W'oo(O) 


^ 3 ( 0 ) 
20-2 ( 0 ) 
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Figure 5. The normalized height see ( |53[ ), as a function of time 
n = 1 ,..., 200 in the case of rational (p^(0) = | (top) and irrational 
ifi'iO) (bottom). The red horizontal lines show the values Q(0) and Q(|). 


is constant and (531 converges a.s. to Q{0). 


Case 2. If <p'(0) 7 ^ 0, then the set of a.s. subsequential limits of (53) is the interval 


{Q{0 ): 0 S [0, i]}. It is exactly this case which is relevant for random trees. 


2.8. Related results, extensions and open questions. Consider the random 
measure 7 r„ on Z counting the particles in the BRW at time n, 


■= ^ Sz. , 
i=l 


Recall that 6z denotes the Dirac delta-measure at z. The moment generating func¬ 
tion of 7r„ is given by 


A„(/3) := [ e^V„(dt)=e‘"(«"IF„(/3), 

Jr 


where Wn is the Biggins martingale defined in (lOl. As n —)■ oo, we have 


(55) 


sup 


e-‘"(«"A„(/3) - IFoo(^) 


0 , 


where C C is a sufficiently small disc of radius e centered at 0. To use the 
terminology introduced recently in m, ED, the sequence of random measures 7 r„ 
converges with probability 1 in the mod-;/? sense. Roughly speaking, (55) says 
that the random measure 7 r„ is close to the n-fold convolution of the measure 
which has cumulant generating function (p{l3)', see ([^. The random analytic 
function Woo describes the “convolution difference” between these two measures 
and therefore plays a central (and somewhat mysterious) role in the theory of mod- 
(/? convergence. Note that in our case, the limiting function Woo is random. Mod-;/? 
convergence is a very strong notion since it implies many classical limit theorems 
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(like the central and the local limit theorem). In fact, Corollary |2.8| could be 
deduced from the general framework recently given in Feray et al. |2i| ( 

Let us stress that in our proofs we do not use the martingale property of Wn- 
We need only that w.p. 1 there is a uniform convergence of Wn to a random ana¬ 
lytic function Woo and that the speed of this convergence is superpolynomial; see 
Lemma |3.3| Results similar to those obtained in the present paper should hold for 
any sequence of “random profiles” : Z —K (not necessarily related to the BRW) 
having a moment generating function that converges (after it has been divided by 
its expectation) superpolynomially to some random analytic function Woo- 

Our results are strong limit theorems in the sense that they refer to the a.s. 
asymptotic behavior of various random quantities related to such as L„(A;„), 
the mode and the height M„. In the case of L„(/c„), it is possible to augment 
our results by more refined asymptotics for the difference between the quantity of 
interest and its limit. To this end, it suffices to take more terms in the asymptotic 
expansion. The a.s. limits of the rescaled differences have non-normal distribution 
that can be expressed through the derivatives of Woo ■ 

Our main motivation was to apply the results to random trees (which we plan 
to do in a separate paper), so we did not seek for optimal moment conditions. The 
moment condition required in Assumption D can certainly be relaxed. It should 
also be possible to extend the results to BRW’s on the d-dimensional lattice Z"^. 
On the other hand, it should be possible to obtain an asymptotic expansion for 
BRW’s satisfying the strong non-lattice condition. Indeed, the first two terms were 
obtained very recently by Gao and Liu [53]. The expansion stated in Corollary 2.8 


Equation (33), is valid uniformly in the range of “central order statistics”, that is 


as long as k € ((p'{l3-)n -j- en, Lp'{j3+)n — en). The case k = '-p'{l3±)n =F o{n) is not 
covered by our results because the function Woo (/3) is not analytic at /3+ and /?_. 


3. Proofs 

3.1. Notation and method of proof. Consider a branching random walk in 
discrete time satisfying Assumptions A-E. Fix some /3 G and recall the notation 

(56) = ^'(/3), a^{l3) = ^"(/3). 

Consider a random measure Vn counting the appropriately weighted particles of the 
branching random walk at time n: 


(57) 


— p-v{P)-r. 


Vn = e 


E' 




Zj^n - Ai(/3)n 
a{^)yfn 


Here we have written (5(z) instead of 5z as this is typographically more convenient. 

Note that Vn depends on /3 but we usually suppress this in our notation. We 
have 

i/„(K) = W„(/3) ^ IFoo(/3), 


where IF„ is the Biggins martingale given by (10). The measure Vn is finite, but, 
in general, it is not a probability measure. 

The characteristic function of the random measure Vn is a random function 
{i/i„(s; /3): s e R} given by 


(58) 


V’n(s;/?) = 


V„(dz) = fn{s]j5)Wn P -b 


IS \ 

T{P)^/n) 
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where 


(59) 


/«(s;/3) = exp<^ -rupiP) - 


o'(^) 


■TK/J /3 



3.2. Expansion of the characteristic function. In Proposition |3.5| we will state 
an asymptotic expansion of '0„(s;/3) in powers of On the formal level, this 

expansion is a product of expansions of /„ and Wn- The main idea is that /„ can 
be represented as a characteristic function of a sum of i.i.d. random variables (so 
that the classical Edgeworth expansion applies), while Wn can be expanded into a 
Taylor series. Multiplying the expansions of /„ and Wn we will obtain an expansion 
of ipn- 

Let us first state an expansion of fn- Define polynomials Pj{z;j3), j G No, by 
the formal identity 


exp 



^J+2(/3) 

(j+2)! 




OO 


'^Pj{z;/3)uP 

j=o 


j has degree 3j and its coefficients can be expressed through 


It is known that 

... ,Kj+2{P)-, see [31 p. 52]. Also 


Pj is an even (resp. odd) function if j is 


even (resp. odd). The first few polynomials Pj are given by 


Po{z-p) = l, Piiz;p) = 




P2{Z\P) = 


Ki{P)z‘^ kI{P)z^ 


24 


72 


h ( l^b{P)z^ K3(^)n4(/3)2:^ kI{P)z^ 

The next lemma is a central element in the proof of the classical Edgeworth expan¬ 
sion; see Lemma 11 in [34l Ch. VI, p. 175]. 


Lemma 3.1. Fix r € Nq and a compact interval K d Pip. 
jsj < and P G K, we have 




3=0 


0'(/3) 


;/3 


- nr/2 


where Si(n) does not depend on s and lim„_,._|_oo (5 i(r) = 0. 


For all n S N, real 

(i+isr+3). 


Proof. Recall that i?i is the intensity of the branching random walk at time n = 1. 
Let Zi,Z 2 ,... be i.i.d. random integer-valued random variables with probability 
distribution 

(60) P[Zi = k] := e-^(^)+'^'=di({/c}), kGl. 

Then the log-characteristic function of Z\ is just 

(‘ 

(61) logEle*®-^!] = (fiiP + is) - ip{P) = ^ Kj{P)PP^, 

i=i 

where the second equality holds for sufficiently small Js]. In particular, EZi = ^(/3), 
VarZi = a‘^)p). Consider the random variable 

_ Zi + ... + Zn — ti{P)n 
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The characteristic function of S* is fn{s]f5); see (59) and (61). The moments of 
Zi are finite because f3 S and is open. Hence we can apply Lemma 11 
from |M1 Ch. VI, p. 175]. In fact, by looking at the proof of this lemma, one can 
see that lim„_>oo Si{n) = 0 holds uniformly in P G K. This yields the statement of 
the lemma. □ 


At this place, we can prove Proposition |2.6 


Proof of Proposition \2.(^ Let Zi,Z 2 ,... be i.i.d. integer-valued random variables 
with probability distribution (60) and log-characteristic function (61). The j-th 
cumulant of Zi is HjiP)- We have 

e'^'=-^(^)"EL„(fc) = = F[Zi + ... + Zr, = k] 

Applying the classical asymptotic expansion ([^ , we obtain (28). By examining the 
proof of the classical expansion, Theorem 13 in [34l Ch. VII, p. 205], one can see 
that all estimates there are uniform in P G K. □ 

The next lemma is a Taylor expansion of Recall that (/3_, /3+) denotes the 
open interval consisting of all P G such that (p'{P)P < ^{P)- 

Lemma 3.2. Fix r S No and a compact interval K C (/3_,/3+). Then we can find 
a sufficiently small £ > 0 such that for all P G K and every u G C with juj < e, 


WooiP- 


W, 


U)/ 


. VVcc {P) j 


1=0 


<MPuY^\ 


where Mi is an a.s. finite random variable. 


Proof. By [9l Corollary 3], for sufficiently small e > 0 the random function Woo(P) 
is analytic on the domain K 2 e ■= {P G C: dist(,5,Ar) < 2e}, with probability 1. 
Fix some P G K. The random function 


r]{u-,P) := 


1 


.r+l 


W^{P- 


^)-E 

1=0 


w^\p) 

j! ' 


0 < Ju] < £, 


can be extended to u = 0 by continuity. It is a.s. analytic, hence the maximum 
principle yields that 


|?7(u;/ 3)| < sup \'q{u]P)\ < [sup \W^{P + u)\ + E 


|«|=e 


|«|=E 


1=0 


The right-hand side can be bounded by an a.s. finite random variable Mi not 
depending on P G K, thus proving the statement. □ 

The next lemma shows that Wn converges to Woo at an exponential rate. The 
proof closely follows the arguments given in Section 3 of [S]; see also [H Section 3]. 
It is here that Assumption D is crucial as it allows us to move into the complex 
plane and to use Cauchy’s integral formula in order to obtain a suitable maximal 
inequality, see (63) below. Such an inequality is the basis for sup-norm convergence 
in the space of continuous functions on a compact set. 
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Lemma 3.3. Fix a compact interval K C {l5-,j3+). Then we ean find a sufficiently 
small e > 0 such that for all n € N, 

sup sup |W'„(/3 + m) - Woo(/3 + m)| < 
where M 2 is some a.s. finite random variable. 

Proof. Let /3 G K. The first part of Lemma 2 in implies that for some p > 0, 
and with 

Zi{t) := f) + 2pe^'^®‘, 0 < f < 1, 

we have that, for some C = C{(3, p) < oo and some k = n^f), p) < 1, 

(62) ^\Wn+i{zi{t)) — Wn{zi{t))\ <Ck^ for all n € N, 0 < t < 1. 

Further, let Dp(/3) = {z G C : \z — fi\< p}. Then, by Lemma 3 in [5], 

(63) sup \Wn+l(z) - Wn{z)\ < f Wn+l{zi{t)) -Wn{zi{t)) 

zGDp(P) JO 


dt. 


Combining (62) and (63) and using Fubini’s theorem we arrive at 

/•i 


E 


sup |VF„+i(z) - VF„(z)| 

zGBp(/3) 


< 


E 


10 


Wn+l{zi(t)) - Wn{zi{t)) 


dt < Cl 


As K is compact we can now find an e > 0 such that 




sup sup |VF„+i(/3 + u) - Wn{l3 + u)\ 

./SGif |ti|<e 


< OO, 


and we may then take 


M 2 ■■= '^e^"'sup sup I W'„+i(/3 + It) - W'„(/3 + m)|. 

/3gK|«|<e 


□ 


n—1 


Using these tools we obtain 

Proposition 3.4. Fix r G No and a compact interval K C (/?_,/?+). Then we 
can find a sufficiently small e > 0 such that for every /? G AT, n G N and s G C 
satisfying |s| < ea-{fi)^/n we have 




1 W^iP ) ( is 

7! 

i=o 


o-(/3)v^ 

where M 3 is some a.s. finite random variable. 


o-(/5) 


< 


M- 


r+ 

n 2 


r(l + kr' 


Proof. Apply Lemmas 3.2 and 3.3 with u = is/{a'{j3)y^n). 


□ 


The next proposition states the asymptotic expansion of the characteristic func¬ 
tion -ipfs-fi). 

Proposition 3.5. Fix r G No and a compact interval K C (/3_,/3+). Then, for 
sufficiently large n G N, /3 G AT, and all real |s| < we have 

S2{n)^_i,2 


(64) Ar{s,n-,l3) := iffs] f) - e Vr,n{s] ft) 


< 


r/2 


e 2 ® (1 + 


|4r+3 


), 
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where 


(65) = “ 


^i/2 jyji 

m=0 


(/ 3 ) 


P, 


j-m 


ZS 

W) 


;/3 


and S 2 {n) is a sequence of random variables (not depending on s or (3) such that 
lim„_>+oo S 2 {n) = 0 a.s. 

Proof. Let e > 0 be sufficiently small to ensure that the statement of Proposition |3.4| 
holds. The subsequent estimates are uniform in (3 G K. We use the inequality 

l^lPl — ^2P2| < \3i.l — ^2||Pl| + |^2||Pl — P 2 I, 

which is an easy consequence of the triangle inequality, with 

is 


Ai = f„is;l3), 

r 

E 


Bi = Wn[ 13 + 


(t{I3)^J ’ 




1 W^\(3) f is ^ 

j=0 




Recall from ( |58| that ifn{s', (3) = ^iPi- Condition |s| < implies that (for sut 
ficiently large n) we have |s| < £a{l3)yjn and hence Lemma 3.1 and Proposition 
yield the estimates 


3.4 




3r+3 


M 3 


), iPi-P2i<^(i + isr+^). 


+1 
n 2 


In particular, since Pj has degree 3j and |s| < we have the estimates 

|P2|<M4, |Pi|<M4, IZI 2 I <Ce-5"', 

where M 4 is a.s. finite random variable and C > 0 is a constant. By the above 
estimates it follows that 


3’'+3)M4 + Ce-5'*'^(l + |sr+i) 

n 2 


(66) |V.„(s;/3)-A2P2|<^e-5«^(l + |s 

where < 53 ( 71 ) is an a.s. finite random variable and lim„_,.+oo <53(n) = 0 a.s. Also, in 
the double sum 


A 2 B 2 - e 


.^(;+m)/2 * I tT(/3) 

1^0 m^O ^ 


IS 


;/3 


wtHP) ( is 


m 


the sum of all terms having j := I + m < r is e 2 ® Pj.„(s;/5). The sum of the 
remaining terms can be estimated by 


(67) 


|A2P2-e-^V,,„(s;/3)| < {1 + Isf^. 

n 2 


The statement of the proposition follows from ( 66 ) and (671. 


□ 
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3.3. Proof of Theorem 2.1 Recall from (18) that 


lik) = 


k — ^i{/3)n 
a{P)y/n 


We can write the characteristic function of see (58), as follows: 




kGl. 


This function is periodic in s with period 27rCT(/3)-y/n. Inverting the Fourier trans¬ 
form we obtain that for fc £ Z, 

1 r+Tr<T{^)^ 

(68) = —J 


^„(s;/3)e-*^""('=)ds. 


' —7ra{0)y/n 

Finally, we are ready to prove the asymptotic expansion of Ln{k). To this end, 
we will insert the expansion of ■)/'„ into (68). Let Z? = ^ be the differentiation 
operator. Define polynomials Qmj, j G No, m = 0,..., j, by 


(69) 


D 


a(/3)’ 


;/3 e 2 ^ = Qm,jiz;/3)e 2 ^ . 


Note that the function qj{z\ f5) appearing in the classical Edgeworth expansion (28) 
is given by 

^ ;/3) 


SO that qj{z;P) = Qoj{z-,P). 

Remark 3.6. It follows that both Qm,j{z]j3) and qj{z;P) are linear combinations 
of Hermite polynomials Hefc(z) with coefficients depending on /3. For even (resp. 
odd) j, all Hermite polynomials involved have even (resp. odd) index k and hence 
the functions Qm,j and qj are even (resp. odd). 

Fix r G No and a compact interval K C (/3_,/3+). Our aim is to prove that 


(71) 
where 

(72) 


sup sup n 2 
peK kez 


a(/3)v^e'5'=-‘"(«"L„(fc)- 




Ur,n (^) ( /?) 


Ur^n{z;P) - y] ,2 ^ 


wt\p) 


77,f/2 rn\ 

j—Q ra—Q 


/5)- 


The idea is that on the right-hand side of (68) we can approximate tjjnis', 13) by 
e“5'’ I4.n(s;/3) and the integration range can be replaced by K. Namely, we will 
show below that 
(73) 

+TTa(fi)^ I- 

V’n(s;/3)e-**^"('=Ms- / e"^* K.n(s;/?)e-*"^"('=M, 

Jm 


n 2 sup sup 
/3Gif feez 


' —■KG(P)y/n 


0 . 


The main tool in the proof of (73) is Proposition 3.5 Given that (73) holds, we 
can replace the integral on the right-hand side of (68) by 

^ f e-2^'K,„(s;^)e-*«""«ds= ^e-2"'WD,,„(x„(A:);/3), 

zT^ Jr v27r 
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where Ur,n{z]P) is as in (72|. To see this, recall (65) and note that for all I G No, 




e 2 ® 




Ms = — 
27r 


e 2 ®^ I — 


D 

W) 


Ms 


(- 


D 


V o-(/3) 


e 2 ^: _ 


Further, 


K/3) 


P 

^3-m 


m 


;/3 e 


is the Fourier transform of the function 


Z i-A 


1 


D 


\ 0’(/3) 


a{py 


/3 e- 


Hence, using ( |69[ ) and Fourier inversion, 
1 


27r 


e 2 ®^ 


K/3) 


p 




Ms = 




Qm,j(z;/3)e 2 2 = . 


In the rest of the proof we establish 
Step 1. It follows from Proposition |3.5| that 


,1/7 


(74) 


sup 


V'n(s;/3) -6 2 ® 14 „(s;/3) 


f3eKj-n^/-r 

Step 2. We prove that for any a > 0 small enough 


ds ^ 0. 

n—>-oo 


(75) 


r/‘I 

n ' sup 


|i/i„(sM)|ds —4 0. 


The functions IF„, n G N, can be bounded by an a.s. finite random variable in some 

sufficient to prove 

l/n(s;/3)|ds ^ 0. 


neighborhood of K. By (58), it is therefore sufficient to prove that 
(76) 


sup 


/SGif Jni/7<|s|<aV?i 

Recall that fn{s',l3) is the characteristic function of {Z^ + ... + Zy)lyjn, where 
Z* := {yZi — p{j3 ))/are i.i.d. random variables with distribution given in (60). 
Note that EZ* = 0, VarZ* = 1, and c := E|Z*p is bounded as long as (3 G K. 
By Lemma 12 from [33] (where we take b = 1/2) we have, for |s| < y/nj{2c), the 
estimate 

sup |/n(s;^)| < e"5®". 

P&K 


This implies (76) with any a < l/(2c). 

Step 3. We prove that for a sufficiently small a > 0, 


(77) 


sup 


|'0n(s;/3)|ds —4 0. 

/3GArdaVS<|s|<7i-<T(/3)VS 
The main difficulty of this step is that, on the range in ( |77| ), we cannot claim 
that the functions Wn-, n G N, are uniformly bounded by some a.s. finite random 
variable. Instead, we shall employ an estimate on the moments of obtained in 
Lemma 2 of |3]. We consider only s G [ayM, 7rtT(/3)yM], the case of negative s being 
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analogous. Write u := s/{a{p)^/n) e [a/cr{^),7r] and z := /3 + iu, so that /3 = Rez. 
By (58) and ( [^ we have 

1/ ( |W"n(z)| 

\W'n\^:P)\ QUifCRe z) — n'Retp{z) ' 

Let ao = inf^gi^(a/(T(/3)) > 0. To establish it suffices to prove that 

(3^K aQ<u<Tz 6 

We shall show that for every zq S C such that Rezo € (/3_,/3+) and Imzo € 
[oq, tt], there are sufficiently small Eoj i^o > 0 and sufficiently large Cq (all quantities 
depending on zq) such that for all n G N, 


j^h+l)/2 gup SUP \Wniz)\ ^ ^ 


(79) 


E 


\W^{z)\ 

.eD^^o) e'^v{Rez)-nRe^(z) 


< Coe 


where ©^^(zo) = {z G C: |z — zo| < Eq}- Indeed, given (791, we can cover the 
compact set K x [ao,7r] C C by finitely many disks of the form Dep(zo), and then 
use Markov’s inequality to obtain (78). 

In the following we prove ( [7^ . We have (/^'(Rezo) Rezo < (/^(Rezo) by the as¬ 
sumption Rezo G {I3-,I3+) and the definition of the interval (/3_,/?+) given in (12), 
(13). Hence, we can find a G (1,2) close to 1 such that oRezo G (/?_,/?+) and 
a(p(Rezo) — (/^(aRezo) > 0. Further, by Assumption E, (p(Rezo) > Re(/3(zo) be¬ 
cause Imzo G [oojTr] and uq > 0. (This is the only place in the proof where we use 
Assumption E). Therefore, we can choose do > 0 such that 

(80) 2do < min{i^(Re zq) — a~^<p(a Re zq), (/^(Re zq) — Re ip{zo)}. 

For z G C such that aRez G (/3_,/3+) define 

_ ^viaRo z)-aRoip{z) 

By the continuity of the functions (p and k, for eg > 0 small enough, the disk 
D3eo('2o) is contained in the strip {z G C:aRez G (/3_,/3+)} and for all z G 
Dseoi^o), 

(81) gTii^(Rez)-nRei^(z) ^ ^nip(Re Zo)-nRe ip{zo)-nSo/3 ^ ^ n(zo) . 

By the Cauchy integral formula, see Lemma 3 in [3], 


(82) 


sup \Wn{z)\ < - 

zGBej,(zo) 


r2TT 


|W„(zo + 2eoe*")|d0. 


Combining estimates (81) and (82) and using Jensen’s inequality afterwards, we 
obtain 


E 


sup 

zGBej,(zo) 


\Wn{z)\ 


0n(/?(Re z) — n Re <p{z) 


< 


/q E|W„(zo + 2goe-^)|dg 

yj-0n(/3(Re zo)—n Re Lp{zo) — nSo/3 

/Q""(E|W„(zo + 2goe-^)|")i/“dd 
Tre 


,n(/?(Re ZQ) — n Re (p{zo) — nSo /3 
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By Lemma 2(ii) from [9] (which uses Assumption D), uniformly over 9 G [0,27r] it 
holds that 


^ n— 1 


1 /a 


{E\Wn{zo + 2eoe*^)!“)'/“ < C' ^ k^{zo + 2eoe^^) 


\r —0 


/ n — 1 


l/o 




\r —0 


where in the second inequality we applied (811. From the above two estimates and 
the definition of k it follows that 


E 


\Wniz)\ 

.gd 3 -o) 


< ^e“ V(aRezo)-v(Re2o) \y gRe v(zo)-¥3(Re zq)^ " < 


where in the last step we used (80). This completes the proof of (79). 
Step 4. We show that 


(83) 


r/2 


f e-5'*>,,„(s;/3)|ds ^ 0. 

7|s|>ni/7 


By (65) we have that for n S N, 

|F,.„(s;/3)| <M(l + |s/’-), 


where M is an a.s. finite random variable. This easily implies (83). 


Taken together, the results of Steps 1-4 imply (73) and thus complete the proof 
of Theorem 12.11 


3.4. Proof of Proposition |2.5| . The idea is to take sums in Theorem |2.1| and 
then to approximate these sums by Riemann integrals. 

Step 0. By Theorem |2.1| with r = 3 we have 
(84) 

L„(h) / Fi(t„(/i);0) F2(t„(/i); 0) \ f 1 \ 

1Too(0)to" (T(0)V27rn V lLoo(0)v^ lFoo(0)n J 

where the o-term is uniform in h G Z and in the formula for x„(h) we take (3 = 0. 

Step 1. By the standard error term analysis for the trapezoidal rule, for any 
integers fci < k 2 and every function / S C^[ki,k 2 ] we have 


k2 


k2 


E/w = + / f{t)dt+- E 


h—k-i 


h—kx-\-l 
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where & [h — 1, h]. Using this formula for the function f{xn{h)) instead of f{h) 
we obtain 


k-2 


V /(x„(/i)) = ^/(a;„(fci)) + i/(a;„(fc 2 )) + cr(0)V« / f{z)(lz 

h=k, 2 2 

12cr2(0)n 

In our applications, /" decays exponentially so that we can estimate the error 
term involving the sum of f"{^h) by 0{l/^/n) uniformly in ki,k 2 - In particular, 
uniformly over |fc| < we have 

(85) 

( 86 ) 

^ rXn(k) 

'Y Fi{xn{h)',Q) = a{0)y/n e“ 2 ^ ^^^( 0 ; o)d 2 : + 0(1) a.s., 

h=-[n3/'t] 

(87) 

k 

Y e“53:^(O^2(a;„(/i);0) = 0(V?T') a.s., 
where we dropped the integrals over the range z < —t?!^ because of 


( 88 ) 


[ e 2 ^ dz = o(n ^), [ e Fi(z;0)dz = o(n a.s. 

7 | z |>[„ 3 / 4 i J\z\>\n^/*] 

Dbtain 

wLio) 


for every O > 0. Recalling the formula for Fi(z; 0), see (21), we obtain 


e 2 *Vi(z; 0 )dz =-^^^^^^^^^^He 2 (a:„(fc))- 


6tT^(0) (t(0) 

Let I A: I < Taking in (84) the sum over h = — ..., fc and using (85), 

(86), (87), ((^, we obtain 


(89) 


IUoo(0)r 


k 1 pXn{k) 

Y e- 5 ^"dz = 

/t=-[n3/4] V “'-2 

1 ^2 ( 




a{Q)^/2'Kn V2 6(T^(0) 


(x„(/c) - 1) - 


W^^(O) 

Woo(O) 


a.s. 


Note that the error term o(n ^/^) in (84) taken n^/'^ times can be estimated by 
o(l/V«). 

Step 2. It remains to estimate X)|ft,|>[n 3 / 4 ] Ln{h). It is an easy consequence of the 
Chernoff inequality that for every O > 0, 


^ ^ EL„(h) = o(n (^+2)). 


|/i|>[n3/4] 
















EDGEWORTH EXPANSIONS FOR PROFILES OF LATTICE BRW’S 


29 


a.s. 


By the Borel-Cantelli lemma and the Markov inequality. 

(90) ^ 

\h\>ln3/^] 

This completes the proof of Proposition |2.5[ 

3.5. Proof of Theorem |2.13[ Recall that, in the first part of the theorem, is 
an integer sequence such that for some a S K, 

= v?^(/3)n + acr(/3)\/n + o(v^), n —)■ oo. 


Note that Xn := Xn{kn) = a + o(l); see p^ . By Corollary 2.7 with r = 1, we have 

(91) = e"^-' • /3) + o f 

y/n Xy'n, 


a.s., 


Here, is given by (31), hence 

Xn 


(92) 


F^XnlP) = 


a 


+0(1) a.s. 


a 


e-5“V;,(/3)+o( ^ 


Inserting (921 into , we obtain 

e^f‘^-‘^(P)-V^a{l3)Ll{kn)= ^ 

This proves the first part of Theorem |2.13| 

Now let kn be an integer sequence such that 

kn = tp'{P)n + Cn, where c„ = o(^/n). 

Note that 


1 

^/n 


a.s.. 


(93) 


Xn •— Xnikn') — 


<x{P)Vn 


= o(l)- 


By Corollary 2.7 with r = 2, we have 
(94) e^'="-‘^(^)"v^o-(/3)L°(/c„ 


-ia;^ 
= e 2 " 


^F°{xn;l3) + -F°{xn;l3) 
Jn n 


a.s., 


where F° and F^ are given by equations (31) and ( |32[ ). Using these equations 
together with (93), we see that the asymptotic expressions for F^ and F^ look as 
follows: 

(95) -^F°{xn;P) = ^^WUl3), 


/n ‘ a^{j3)n 


(%) = 


WM- 


WM]+o[-] a.s. 


2aHP) 

Here we used that He4(0) = 3 and He2(0) = —1; see ([^, ([^. Also, note that 
by (@, 


(97) 

Now we consider two cases. 


e =1 + o(l). 


Case 1. Let c„ —>■ +oo but c„ = o{^/n). Inserting (95), (96), (97) into (94), we 
obtain 

e^'="-^(«"y^cT(/3)L“(fc„) = ^^lU^(/3)+o(^) a.s. 

f7^(p)n \ n / 
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This proves the second part Theorem |2.13| 


Case 2. Let c„ = 0(1). Inserting ( pl^ , (96l, (97) into (94), we obtain 
(98) e^'="-‘^(«"v^a(/3)L“(A:„) 


«^3(/?) \ 

n Wcr'^W) ' 2(7^iP)) 


+ 


wM - 


2a2(/?) 




01 - 


a.s. 


This proves the third part of Theorem |2.13[ 


3.6. Proof of Theorem 2.17, We start by proving the formula for the mode 
Un stated in the first part of the theorem. We search for k G Z maximizing the 
occupation number Ln{k). 

Step 0. We write k = kn{a) = ^p'{Q)n + a, where a is subject to the restriction 
a G Z — ip'{0)n. Let /3 = 0 and = (p''{0). With this notation, we have 

Xn ■= Xn{k) = — 


By Theorem |2.1| with r = 2, we have 
(99) 


m 


V27rn tTL„(fc) = e ( Woo(O) + ^Fi(a;„; 0) +-T' 2 (a;„; 0))+o 
\ i/n n J 


a.s., 


where Fi and F 2 are given by equations (21) and (22). The o-term is uniform over 
a G Z — ip'{0)n. 

Step 1. Let us first assume that |a| < for some small e > 0. Later we 

will show that, asymptotically, the values of a outside this range have no chance 
to be the mode. Let us agree that all estimates will be uniform over |a| < 
Recalling (21) and (22) and noting that = o(n“^/^) and x„ = o(l) we obtain 


( 100 ) 

( 101 ) 

where 

( 102 ) 


1 


«3(0), 


—Fi(x„; 0 ) = - _W^( 0 )-^Woo( 0 ) +0 - 


/n n \ a 

-F2(x„;0) = -C + 0 
n n 


2cr4 


a.s.. 


a.s.. 


C = 


^ 4 ( 0 ) 5 ^ 3 ( 0 ) 


24cr6 


W^oo(O) + 


and we used that He4(0) = 3, He6(0) = —15, He2(0) = —1; see 0, 0. Noting 
that = o(l/n), we obtain by Taylor’s formula that 


(103) 


e =1 — 


2cr^r 


01 - 


Inserting (100), (101), (103) into (99), we obtain 
(104) 


m crL„(fc) 

= W^oo(0) +Va( ^W^(O)- 


^ 3 ( 0 ), 

2cr4 


W^oo(O)) -^IRoo(0) + C) + 0 


a.s. 
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Differentiation shows that the arg max of the quadratic function in the brackets is 
given by 

wLio) « 3 (o) 


(105) 


= 


W'oo(O) 2 ct2 ■ 


However, we have the restriction a G Z — (p'{0)n. Using (104) we easily obtain that 


(106) 


l'ima{Ln{k + 1 ) - Ln{k)) = 


- -- a 


01 - 


It follows that there is an a.s. finite random variable such that for n > Ni, the 
expression in (106) is negative for a > a* — | and positive for a < a* — Hence, 
for n > Ni, the mode (computed over the range |a| < is equal to one of 

the numbers [u*J or where u* = (p'(0)n + a*. 

However, note that we do not claim that the mode is equal to w* rounded to 
the nearest integer! Indeed, if m* is so close to a half-integer that the first term on 
the right-hand side of (106) is smaller than the error term o(I/n), then we cannot 
tell whether [u*J or [u*] is the mode without considering further terms in the 
asymptotic expansion. 


Step 2. Inserting into (104) any a = 0(1) we obtain that 
(107) ^ .. 


m "^27171 (jL„(fc) = Woo(O) -I- O 


a.s. 


Take any B > 0. In order to complete the proof of Theorem 2.17| (a) it suffices to 
show that 


(108) 


B 


lUoo(O) — m ^V2Trn aLn{k) > — 


for all n > N 2 and |a| > . Here, we denote by iV 2 ,iV 3 ,... a.s. finite random 

variables. The proof of (108) will be given in Steps 3-5 below. 


Step 3. Let |a| > a^/n. Then |a;„| > 1. By Theorem 2.1 with r = 0 we have 
(109) m“"-\/27rncrL„(A:) = e“5'^"IUoo(0)-I- 0 ( 1 ) a.s. 

The o-term is uniform, hence 

IToo(O) - m~^V2TTnaLn{k) = IToo(0)(l - e“5“") - o(l) > ^IUoo(O) > — 

6 n 

for n > N 3 , where we used the inequality 1 —e“5 > 1 and the fact that lUoo(O) > 0 

a.s. 

Step 4. Let |a| < ay/n but |a| > for some e > 0. This means that |a;„| < 1 

but \xn\ > By Theorem 2.1 with r = 1 we have 


(110) m '^V2TTnaLnik) = e 2 “" ^IUoo(O)-|-^Fi(x„; 0)^ -I -1 


1 


a.s. 


Recalling (21) and using the condition |a;„| < 1, we obtain that T’i(a;„;0) = 0(1). 
Hence 

r.2 


IUoo(0)-m-”v^aL„(A:) = IUoo(0)(l-e-5"')-0 > '^W^{0)-0 


1 
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where we used the inequality 1 — e > yji valid for y S [0,1]. By the inequality 
the right-hand side is bigger than B/n for n > N 4 . 

Step 5. Finally, let |a| < but |a| > . Equivalently, |x„| < 

but \xn\ > By Theorem 

( 111 ) 


2.1 


with r — 2 we have 


V27rn aLn{k) = e 2 ®" (iFoo(O) + ^Fi(x„;0) -I- -F'2(x„;0) ) -I- 
\ y'n n J 

The above assumptions on x„ together with (21) and (22) imply that 


01 - 
n 


a.s. 


1 


1 


^Ei(x„;0) = 0(n _i^2(x„;0) = O - 


/n 


It follows that for sufficiently small e > 0, say e < 1/20, 

2 

VFoo(0)-m-"v^aL„(A:) = IFoo(0)(l-e-5"' )-0 > ^W^{0)-o (xl) , 

where we again used the inequality 1 — e“^/^ > y/3 valid for y G [0,1]. By the 
inequality x^ > , the right-hand side is bigger than B/n for n > N^. This 

completes the proof of (108). 

Step 6 . Finally, we prove Theorem 2.17| (b). We are interested in the height 
Mn = Ln{un), where Un is the mode equal either to [m*J or to [u*]. The idea is 
to insert both numbers into expansion (104). Using ( |104 ) with a = a* + 0, |0| < 1, 
and recalling that a* is given by (105), we obtain 

a Ln{u*^ + 9) IFoo(O) 


- lUoo(O) -h (2(a* -h 0 ) 0 * - (a* + 9^) + — + o 


= VFoo(0) + ^^(a?-fl2) + -+o(- 


C 

n 


a.s. 


Recalling the formula for C, see (102), and after some straightforward transforma¬ 
tions, we obtain 


27rucrL„(u* -I- 0) 


= 1 - 


1 


2cr2r 


^ 3 ( 0 ) K4(0) 


6 a'^ 


4cr^ 


+ 0^ + (log Woo)^^(O) 1+0 ( — ) a.s. 


Now, the mode Un is equal to u* -I- 0 with 0 being either [u* J — u* or [u*] — u* 
It follows that 


1 ^ 00 ( 0 ) 771 ’ 

7 , the m 
)llows tl 
y/27rna Ln(un) 
lUoo(O) to” 

where 9„ is the smaller one of the numbers u* — [u*J or [u*] — u*. 


= 1 - 


1 


2 a^T 


K§(0) K4(0) 


6(7^ 


4cr^ 


+ + (log Woo)"(0) ) +0 ( “ ) a.s., 
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